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A  aodel  tacrlbiig  the  Internal  i>|ii4  cs  of  a  aoMModated  llqnld  has  been 
ooMtivctci  Ter  mc  la  the  eel relation  of  cold  ewtna  iratterlag  cross 
senicBs  mi  free- radical  reccsfclaatlaa  rates  la  liqald  iijfcogu.  The  propose J 
aodel  as* aars  that  the  sped fl cat loa  of  a  glsea  tenperstnre  aad  density  ladles 
a  qvlescest  state  ia  shldb  each  aolecale  vibrates  with  nsdos  phase  Is  a 
field  set  ap  by  the  ia  of  the  poaltlaa  probability  deaaltles  of  its  Belabors , 
the  oeatral  sites  of  shldb  e  nap  rise  a  cloae  pached  ^asal-latt&oe.  A  ooapater 
progras  has  bees  srittea  to  calealate  the  eoaslsteat  vibration  freqweacy,  the 
aolecalar  sell  depth,  the  lstereal  tsergf,  the  eat  ropy .  aad  the  free  energy 
of  a  Hqsld  with  a  gives  1 1  ape  rat  i«s  aad  density  asd  with  tesoasi  Lesnard-Jcoes 
paraarters  of  the  late  molecular  lateraetloa.  Application  of  the  prograa  tc 
the  case  of  ilqsid  syingn  at  points  an  the  P-f-T  surface  rsagH  if  fros  2C* 
to  2ft  alrag  the  sataratlaa  llae  has  glss  predict  lass  of  themodyaaslc 
properties  which  are  la  good  agnoast  with  esperiarstal  data.  It  is  tentatively 
oaacladed  that  the  aodel  anticipates  the  aoat  si^dflemt  geseral  fp stares 
of  sotioB  In  the  llqnid  sad  aqr  be  laplsaastcd  la  ealralatlacs  of  sram* 
sections  sad  recaablaatlas  rates  which  are  depea  dent  an  the  details  of  tae 
late  real  ftmnsl  cs  and  adcroscnplc  strartw*. 
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In  the  design  of  hydrogen-prope  1  led ,  noeleer-pCMersf 
rocket  aytteas,  two  independent  problems  ere  encountered 
which  require  a  mathematical  description  of  the  dynamics 
and  distribution  cf  molecules  in  liquid  hydrogen  for  their 
solution*  First,  it  has  been  shown  (kef.  1)  that  the  molecular 
species  produced  by  the  attenuation  of  neutrons  and  gams 
rays  in  liquid  hydrogen  say  not  immediately  undergo  the  exotLier- 
ral  reactions  leading  to  a  return  to  the  1^  ground  state. 

Su  potential  energy  that  ia  stored  temporarily  in  the  formation 
of  radiation  products  such  as  H,  and  (vibrational  ly 
or  rotationally  excited  molecules)  nay  not  appear  as  heat 
for  several  seconds,  even  assuming  that  the  nobility  of  such 
products  in  the  liquid  is  ss  high  as  in  a  gas  of  the  same 
density.  Since  a  large  fraction  of  the  radiation  energy 
goes  initially  intc  the  formation  of  such  species  and  since 
the  most  intensely  irradiated  portion  of  the  liquid  propellant 
is  generally  nearest  the  exit  of  the  storage  tank,  it  is 
possible  that  the  delay  in  equi partition  will  effectively 
reduce  the  amount  of  heat  production  while  the  rocket  ia 
under  power. 

*?hen  the  lattice  structure  of  liquid  hydrogen  is  taken 
into  account,  it  nay  be  found  that  the  inhibition  of  nobility 
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associated  with  aolecalar  ordering  and  cryogenic  toperatore 
results  in  a  delay  in  the  equi parti tion  of  energy  that  is 
roch  longer  than  seconds.  In  the  case  of  one  particular 
kind  cf  rotationally  excited  molecule,  namely,  ortho-hydrogen, 
the  ret  cm  tc  the  para  -  hydrogen  ground  state  is  known  to 
require  months.  The  evolution  of  heat  over  a  long  period 
presents  a  problem  in  space  vehicles  that  east  retain  part 
of  the  irradiated  propellant  for  future  operations. 

A  second  problem  requiring  the  use  of  a  liquid-hydrogen 
model  is  the  production  cf  hydrogen-capture  gasna  rays  by 
cold-neutron  absorption.  Owing  to  interference  effects, 
the  neutron  scattering  cross  section  of  the  free  hydrogen 
molecule  falls  off  drastically  below  3.023  ev  (Be?.  2). 

Under  the  assurption  that  the  cross  section  of  molecules 
in  the  liquid  is  the  sac e  as  for  free  molecules,  the  cold- 
neutron  man  free  path  in  the  liquid  be  cones  comparable  to 
the  fast-neutron  penetration  distance,  and  the  probability 
of  leakage  fren  the  propellant  tank  ray  become  greater  than 
the  probability  of  neutron  capture.  This  establishes  the 
importance  of  cold-neutron  transport,  since  hydrogen-capture 
gammas  are  inport  ant  fren  botn  the  heating  and  the  dose-rate 
standpoints.  It  follows  that  the  influence  of  liquid  forces 
on  the  oold-neutron  scattering  cross  section,  an  effect  wnich 
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is  sure  to  be  significant  but  is  at  present  of  unknown  nagnitude, 
nas  an  important  bearing  on  shield  and  pressure-vessel  design. 

The  treatment  of  the  two  p rob leas  cited  above  requires 
a  mathematical  description  of  both  the  neighbor  distributions 
for  a  typical  molecule  and  the  state  of  notion  (or  possible 
states  of  notion)  of  the  molecule  itself,  'fith  this  information 
the  nagnitude  of  the  energy  barriers  which  inhibit  the  free 
migration  of  radiation- induced  molecular  species  can  be  determined 
and  the  effect  of  liquid  binding  forces  and  lattice  structure 
on  the  molecular  cross  section  can  be  calculated. 

The  present  report  describes  a  model  for  the  internal 
dynamics  and  microscopic  structure  of  a  nonassociated  liquid 
which  will  lend  itself  to  calculations  of  the  types  required 
in  the  consideration  of  an  irradiated  hydrogen  propellant. 

Evidence  for  the  validity  of  the  model  is  given  in  the  fcm 
of  a  comparison  between  predicted  thermodynamic  properties, 
calculated  with  the  aid  of  an  ISM  7093  computer,  and  published 
experimental  results.  Applications  of  the  model  to  the  problem 
of  cold-neutron  scattering  and  recombination  rates  are  the 
subjects  of  another  NARF  report  (Ref.  3) . 


3 


2.  THE  PROPOSED  LIQUID  MODEL 


2.1  Description  and  Assumptions 

Any  theory  of  the  liquid  state  with  a  claim  to  rigor 
must  be  derived  from  a  general  statistical  mechanics  formulation 
and  must  be  amenable  to  the  inclusion  of  quantum  mechanical 
effects.  In  such  a  first-principle  approach,  exemplified 
in  the  nonquantum-mechanical  case  by  the  work  of  Kirkwood 
(Ref.  4) ,  the  structure  of  the  liquid  is  determined  from 
the  free  interactions  of  many  molecules  and  is  not  assumed 
to  be  locally  similar  to  tnat  of  a  corresponding  solid. 

However,  the  integral  equations  which  result  are  difficult 
to  solve  without  the  introduction  of  mathematical  approximations 
that  tend  to  offset  the  rigor  of  the  physical  foundation. 

Hence,  nonrigorous  models  based  on  a  prior  assumption  of 
a  •quasi-crystalline"  liquid  structure  continue  to  be  used 
because  of  their  heuristic  value  in  des.  ibing  the  general 
features  of  internal  liquid  motion  and  because  of  their  adapta¬ 
bility  to  the  correlation  of  experimental  data.  It  is  a 
model  of  thi3  sort  that  is  introduced  here  as  a  basis  for 
specialized  types  of  calculations. 

In  the  theory  of  Lennard-Jones  and  Devonshire  (Refs.  5 
and  6)  each  molecule  in  a  liquid  is  assumed  to  move  in  the 
intermolecular  potential  field  that  would  be  set  up  if  its 
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nearest  neighbors  were  *rmeared  uniformly  over  a  spnerical 
surface  passing  through  their  central  lattice  sites.  An 
extended  form  of  the  L-J-D  model  takes  vacancies  into  account 
by  minimizing  the  Helmholtz  free  energy.  However,  the  ratios 
of  vacancies  to  molecules  predicted  in  this  way  are  found 
to  be  much  smaller  than  those  obtained  from  X-ray  determinations 
of  the  liquid  structure  (Refs.  7  and  8).  It  is  inferred 
that  the  latter  model  cannot  be  relied  upon  to  provide  a 
liquid  description  adequate  for  the  calculations  of  quantities 
which  are  dependent  on  the  details  of  the  internal  dynamics 
and  nicroscipic  structure.  The  sane  conclusion  is  reached 
in  regard  to  the  currently  eminent  "significant  structure" 
theory  (Ref.  9) ,  which  provides  good  thermodynamic  results 
but  which  relies  on  the  use  of  an  empirical  Debye  temperature 
and  prescribes  no  method  of  determining  the  potential  energy 
of  the  quasi-lattice  *t  points  off  the  vapor-point  line. 

The  object  of  the  present  effort  is  to  adumbrate  a  tractable 
liquid  node!  which  is  core  or  less  in  the  spirit  of  the  L-J-D 
model  but  which  takes  into  account  both  neighbor  displacerents 
and  quantum  rechanical  effects.  As  in  the  L-J-D  model,  a 
specific  quasi-crystalline  structure  with  local  order  extending 
out  to  the  third  or  fourth  nearest  neighbor  is  assumed. 

However,  in  the  method  offered  here  the  motion  of  each  molecule 
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is  determined  by  the  state  of  action  of  its  neighbor,  so 
that  factors  such  as  the  change  in  molecular  environment 
with  temperature  variation  at  constant  density  are  taken 
into  account. 

The  particulars  cf  the  present  rodel  will  be  formulated 
on  the  basis  cf  the  following  set  cf  as sorptions: 


1.  The  liquid  structure  is  quasi-crystalline, 
with  local  order  of  central  sites  extending 
for  at  least  several  ssolecular  separations. 

2.  At  constant  temperature  and  density,  increases 
in  the  vacancy  fraction  (hole  concentration) 
result  in  a  scale  compression  of  the  lattice 
dissensions  so  that  the  fon  of  the  site 
distribution  is  preserved  with  none  of  the 
molecules  occupying  off-site  positions. 

3.  The  interactions  are  due  tc  central,  isotropic, 
two-molecule  forces  (thus,  the  rodel  is 
limited  to  'nonassociated*  liquids) . 

4.  Each  molecule  roves  quantun-nechanically 
about  its  own  central  site:  ir.  the  vicinity 
of  this  site  the  molecule  roves  "harmonically" 
in  discrete  energy  states  characteristic 

of  a  three-dimensional  harmonic  oscillator. 

5.  Each  molecule  vibrates  in  a  field  set  up  by 
a  "molecular  cloud"  defined  by  the  sum  cf 
the  position  probability  densities  of  its 
neighbors  that  are  averaged  over  ail  modes 
of  vibration;  tne  "consistent"  frequency 

thus  defined  is  common  to  all  cf  tne  molecules, 
as  in  the  Einstein  model,  but  is,  in  this 
case,  determined  uniquely  by  the  common 
neighbor  environment. 

6.  The  foregcing  assurptions  define  a  quiescent 
liquid  in  which  there  is  no  correlation 
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between  the  phases  of  viaratioa  of  neigh¬ 
boring  molecules. 

7.  The  quiescent  systetr  thus  defined  is  disturbed 
by  longitudinal  acoustical  vases  with  a 
velocity  determined  by  the  "consistent* 
frequency  of  vibration  and  with  wavelengt  in 
restricted  to  integral  nultiples  of  tvice 

the  lattice  dimensions. 

8.  Owing  to  absorption,  o.*.ly  those  acoustical 
vases  are  possible  whose  energies  are  less 
than  that  required  for  the  excitation  of 
"consistent"  guantv  levels  is  the  quiescent 
liquid;  the  "ground  state"  energies  of  all 
sash  waves  taken  together  are  identified 
with  those  of  the  stationary,  consistently 
vibrating  nodes. 

9.  For  constant  temperature  and  density,  the 
hole  concentration  tends  toward  a  value 
corresponding  to  the  lowest  Helrholtz  free 
energy;  the  hole  ccncentratioc  is  the  only 
significant  independent  variable  in  the 
ainiuzatioo  of  the  Seinbolti  free  energy. 


Each  of  the  above  assurptions  can  be  supported  by  arguaents 
based  either  on  empirical  considerations  or  on  tbernodynazrical 
verities.  However,  the  ultimate  argurent  for  such  a  construct 
rests  in  its  pragratic  validity,  i.e* ,  its  general  utility 
in  the  satisfactory  prediction  of  empirical  data.  Hence, 
the  above  assurptions  will  be  incorporated  into  wa these t i cal 
statements  without  individual  analysis,  and  the  results  of 
a  subsequent  application  to  the  case  of  liquid  hydrogen  will, 
for  the  present,  be  taken  as  an  indication  of  their  collective 


validity. 


2.2  Coesisteat  Frequency  a^d  Sell  Depth  of  the  Cuies cent 

Liquid 

According  to  Asempticn  1,  the  neighbor  situs  of  each 
aoiccalt  lie  os  concentric  spherical  surfaces  (shells)  ,  with 
the  Mebcr  of  sites  per  shell  determined  by  the  type  of  quasi¬ 
lattice.  Let  and  be  the  radius  and  nuaber  of  neighboring 

nolecules,  respectively ,  corresponding  to  the  i***  shell. 

Let  the  density  of  the  liquid  be  a  and  let  the  hole  concentra¬ 
tion,  defined  as  the  fraction  cf  sites  which  are  unoccupied, 
he  f.  Then,  according  to  Assurption  2,  is  proport icmal 
to  (a/ (1-f )  1  and  is  proportional  to  (1-f) .  The  neighbor- 

site  distributions  out  to  the  fourth  neighbor  shell  have 
been  derived  for  four  types  of  lattice  structure  and  are 
given  in  Table  1.  The  neighbor-shell  radii  conputec  frur 
the  table  are  in  cr  when  M  is  the  molecular  mass  in  aru  and 
:.'A  is  Avogadro *  s  nuaber. 

As  a  setter  of  mathematical  convenience,  the  approximation 
is  nade  (as  in  the  L-J-0  method)  that  the  central  neighbor 
sites  are  "smeared*  over  the  neighbor  shelly  thus  the  environ¬ 
ment  of  each  molecule  is  forced  by  a  radially  symmetric  cloud 
nade  up  of  contributions  from  the  stationary  wave  functions 
of  various  neighbor  molecules  whose  central  siies  are  smeared 
over  the  shells  to  which  they  are  attached.  Figure  1  shows 
the  geometry  associated  w*th  a  central  site  at  0  and  a  neighbor 
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Figure  1  Geometry  Aar.ocialed  with  a  Central  .'tile  a*. 
0  ami  a  Neighbor  .Shell  of  Kailua  a. 


shell  of  radios  r .  The  associated  pocitioc  probability 

dessitv  at  a  distance  r'  fror  9  is  desisted  ms  z  (r'l. 

i 

Then  for  a  spall  displacement  **  of  the  central  molecule, 
the  potentxal  energy  doe  to  aolecnles  in  a  shell  of  thickness 
dr'  associated  with  the  i^  neighbor  shell  is 

•  Ma'.r'Jdr'  *  zi  (r')dr'j  v(s) 2*r'2  sinad#  (2-11 

'l-o 


where 


s 


/r2*a  '2- 


2r's 'cosy 


(2-2» 


and  -rtere  »{s)  in  the  in  te  molecular  interaction.  In  scat 
cases  of  interest  ▼{»)  is  excreuil'ic  in  tae  Leanard-Jcnes 
form,  viz., 

«12 

vCz’  -  -*  r  r  **  (2-31 

c  c 


Carrying  cot  the  integration,  expanding  the  result  xn  terns 
of  o',  and  discarding  ten  of  order  higher  than  s'*  (by 
Assumption  4)  gives 


tI-4) 


Differentiation  with  respect  to  s'  gives  the  restoring  force 
per  unit  central  aolecule  per  unit  displacement  due  to  the 
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itial  shell 


r*  cad  dr' : 


ti*r  >‘*r'  *  [  ]/' 

-  •*r-Vcr'5dr'[5»cr'“#^22»cr-*14  ] 


(2-5) 


Finally,  the  "restoring  coefficient"  doe  to  all  of  the  neighbors 
around  a  central  molecule  is 


K 

1 


(2-*) 


All  that  is  needed  to  evaluate  the  restoring  coefficient. 


and  hence  the  frequency  imposed  on  a  central  nolecule  by 
its  neighbors  in  their  own  quantum  mechanical  distributions, 
is  an  express! ca  for  x^Cr'),  the  radial  density  associated 
niti  neighbors  whose  central  sites  lie  on  the  i1*  cell. 

This  is  obtained  by  integrating  the  density  contributions 
free  all  neighbors  with  central  sites  aaneared  on  the  i 
shell  ac  an  arbitrary  point  is  the  vicinity  of  the  shell. 
Such  am  integration  requires  the  knowledge  of  the  position 
probability  density  of  neighbors  about  their  central  sites, 
which  in  turn  depends  oo  the  "restoring  coefficient"  (or 
classical  frequency)  of  the  neighbors.  Let  (s )  represent 
the  position  probability  density,  averaged  over  all  scodes 
of  vibration,  of  a  neighbor  aolecule  at  a  distance  9  (Fig.  1) 
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from  its  central  site.  (;  will  sot  be  oocfucd  with  t 


liquid  density,  designated  by  Lb 
tc  s  is  indicated  bat  is  aeitted 
ecoal  to  tbe 


for 


**  i  wsbscript 
.)  »^U)  is 


to 


3 

energy  state  weighted  by  tbe  papolatian  fraction  of  socb 
states.  Solatia**  of  tbe  Scbrodiager  equation  with  a  three 
dineesional  hinonic  potential  gives  for  tbe  normalized  wave 
fuactico  of  a  state  characterized  by  tbe  qnantm  stnrbers  (, 
and  n 


C  (z.y.z)  -  s  »  s  H  U*)H  (,y)H  (.ale'172*  *  (2-?) 

i,*,n  leal  **  n 

where  separation  iatc  the  Seirite  polynomials  9  has  easily 

II 

been  effected  by  the  use  of  Cartesian  coordinates  defined  by 


*2.t2.z2 


(2-S) 


The  parameter  a  is  related  to  tbe  stiffness  coefficient  £ 

s 

of  tbe  neighbor  molecnles  by 


•i2.2/)* 


(i-»> 


where  ft  is  the  nolecalar  rass  and  k  is  tbe  Boltzmann  constant; 
the  normalization  factors  are  given  by 


S 

n 


(2-13) 


13 


For  a  Molecule  is  the  2^  energy  state 

l*r*o  *  ^ 

and  tht  nornalized  position  probability  density  is 


(2-11) 


(2-12) 


g^  -  (1/2)  ( j*l)  (j*2) 


(2-13) 


is  the  degeneracy  (cr  nunber  of  terns  in  the  smutioc}  for 

a  three  dinensiecal  hamonic  oscillator.  Then,  if  is 

ti* 

the  fraction  of  nolecoles  m  the  j  energy  state,  the  average 
position  probability  density  xs 


is  5 


(2-14) 


where  for  a  temperature  T 


r 


* 

J 


-E7  kT 

(j*l)(j+2)e  * 


The  energy  levels  E^  are  giver,  by 


E. 

3 


ke 


(2-15) 


(2-16) 
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where  0  is  uniquely  related  to  a  by 


e  «  hV/xk 


(2-17) 


quantity  9  is  related  to  the  classical  frequency  w  by 


hv 


(2-18> 


If  all  of  the  nolecules  vibrate  with  frequency  v,  then  0 


say  be  described  as  the  characteristic  (or  "Einstein") 


perature  of  the  system.  In  tens  of  9, 


-  (j*l)(j+2)e 


'7? 


(j+l)(j*2)e 


-j0/T 


k  (j^l)(j^2)e"j8/7|l-e“#/T) 


(2-19) 


Let 


-e/T 


(2-20) 


Then 


j  3 

■  ia  (1“«) 


and 


(2-21) 


2 

*  (c) 


2  2  2  3 

a  -a  p  --  ■* 


,3/2 


-  e  -  "  (l-«) 


i  V  n» 


j  l,n,r» 


2  4 1 1 


^  (ay)  an  («) 


n  n 

2  n!  2  nl 


(2-22) 
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If  a  and  n  are  fixed  in  the  second  summation  so  that  l  *  j-ra-n, 
then 


♦2(») 


•3  -V  3 

-37J  e  (1-.) 


»V<.y>  H„2 (0*» 


V  ”» 

2nnl 


n 

2  nl 


m 

■z 

j*n+n 


w 


3  H  i-tt-n1"*1 
j-n-n 

2  ( j-ra-n) 1 


•  -a2.2..  .  3  V’  *m  <•*>  Hn  <“>  V 

_  e  <1— »  >  — - V—  > 

■  T  2 2  nl  1=0 


t+m+n  Hl 


2%! 


a3  2  2  t— 1 1  ~  2.  ..2 


3/2 


-•  p  (l-w)3  ura+n  ^  Hn 


ra,n 


m  n 

2  ml  2  nl 


es 

z 

l*o 


.  *  h4  <**> 
III  %  I 


(2-24) 


Application  of  Mahler' s  formula  on  Hermite  polynomials  (Ref.  10) 
gives 

3  2  2 
0  a  -a  0 

♦(p)"^72e  (1'“ 

W 

H  2  (ay)  H  2  (ez) 
n 


|.  2|*1/2 

’.  2  2  -  2  2  2  " 
2a  x  w-2a  x  w 

|l-ui  1  exp 

!-w 

Ex  H  ( 
um+n  n 


n 


m,n 


2  ml  2  nl 


(Equation  continued  on  next  page) 
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a  2(mz) 


,3  22  1 1  f.  «  -  jit-*  a  _c« 

8  *1-->3U--2)  -n 


3/2 


XX 


K  2{«y) 

IE 

rl 


7**0 


a 

"V2 


2tf2 

-a  p 


( !-• ) 3 ( 1-m 2 j  exp  £  j  2  2 ( x2*y2 ) j 


n  E  («z) 


Z»  n  c 

Ifl  -*a 


n*o 


(2-25) 


Thus,  the  position  probability  density  of  a  neigobor  molecule 
with  a  classical  frequency  defined  by  a,  averaged  all  nodes 
of  vibrations,  is  simply 


*2(p) 


-3/2 

» 


3 

X 


.2  2 
-x  p 


(2-26) 


where 


2 

X 


l 


2 

a 


(2-27) 


The  radial  density  distribution  around  a  central  molecule 
can  now  be  calculated  from  Equation  2-26*  Reference  to  Figure  1 
shows  that  the  density  due  to  the  itl1  shell  at  a  point  P  which 
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is  located  a  distance  r '  fro r.  0  is 


x.(r')  *  - 9.2  U(f,r')J  2* a. 

1  4..  2  1 

i  J 


sin  Idl 


where 


•  i 


U %e»r') ]  sin 


(2-24) 


9  *  /a.  ♦r '  -2a  r '  cos  § 
»  i  i 


3y  setting 

.  _  2 


9 

,  i 

vH 

r 

l2  M  2v'  l 

i 

f 

la  -r'l 
1  i  1 

(2-29) 


(2-30) 


-  -[i  *i(  ir^)‘  ^h2(Vrf]  *«•  [-‘2|vrf]j  (2': 

it  can  easily  be  shown  that  the  first  itea  in  Equation  2-31 
is  negligible  for  all  cases  of  interest.  Eance 


x^fr')  -  * 


-3/2  ‘i 


4a.r'  1  *  C3tP 


[-*2lvr1*] 


(2-32) 


The  combination  of  Equation  2-5,  Equation  2-32,  and  Equation  2-6 
gives  for  the  "stiffness  coefficient"  of  a  central  molecule 


(2-33) 


{2-34 


krlTi)  •  l2*c  ( 

and  where  a  -  lr'.  L{t ^ j  is  s  suitably  chosen  lower  licit 
which  oust  be  greater  than  zero  owing  to  the  error  in  the 
neighbor  wave  functions*  assas«d  to  be  hamooic*  at  points 
distant  froa  the  central  neighbor  sites. 

The  well  depth  of  the  central  no  locale*  vf  is  calculated 
in  analogy  to  KR,  using  Equation  2-4  with  o'  set  equal  to 
zero.  Multiplying  the  result  by  Avogadro's  nunber  gives 
the  "nolar  well  depth*"  which  turns  out  to  be 


(2-37) 


(2-38) 
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anc 


The  integrals  in  Equation  2-35,  Equation  2-36,  Equation  2-33, 
and  Equation  2-39  become  infinite  if  L|y,J  is  chosen  to  be  zero, 
this  is  due  to  the  fact  that  the  neighbor  wave  functions, 
vriich  in  keeping  wit.n  Assumption  4  are  harrcnic  near  the 
neighbor  shell,  are  finite  even  at  the  site  of  the  central 
molecule.  The  actual  neighbor  wave  functions,  of  course, 
rust  vanish  at  this  point.  In  practice,  as  long  as  Assurption  4 
holds,  this  effect  creates  no  difficulty.  It  is  found  that 

all  cf  the  significant  contributions  to  k  ,  '<  ,  u  ,  and  u 

a  r  a  r 

core  fro r  portions  cf  the  neighbor  distributions  which  are 
relatively  close  to  the  neighbor  shell  and  far  from  the  central 
site.  In  all  realistic  cases  the  contributions  decrease 
rapidly  with  decreasing  u(*Xr') ,  with  a  large  portion  of 
the  space  inside  the  neighbor  shell  contributing  a  negligible 
amount  tc  the  integral,  f  minimum  is  eventually  reacned 
after  which  the  contributions  increase  drastically  with  decreasing 
r'.  Assurption  4  implies  that  the  harmonic  wave  functions 
are  accurate  down  to  values  of  r'  coresponding  to  negligible 
contributions  to  the  integrals.  Hence,  for  the  sake  of  definite¬ 
ness,  l|y;|  can  be  defined  as  the  largest  cf  the  four  values 
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of  u  corresponding  to  the  minima  of  the  integrands  in  Equation  2-35 r 
Equation  2-36 ,  Equation  2-38 ,  and  Equation  2-39.  The  lower  limit 
is  thus  defined  as 

1ti)  *i[Ti  -/?-"]  (2-40> 

The  salient  feature  of  the  present  model  is  the  require¬ 
ment  that  Assumption  5  be  satisfied;  that  isf  the  frequency 
at  which  all  the  neighbors  of  a  central  molecule  vibrate 
must  correspond  to  a  density  distribution  which  imposes  that 
same  frequency  on  the  central  molecule.  Let  the  "stiffness 
coefficient*  implicit  in  a  frequency  of  neighbor  vibration  v 
(uniquely  related  to  a  through  Equation  2-17  and  Equation  2-18) 
be 

K  -  h2a4/l  (2-40A) 

N 

Then  the  state  of  the  quiescent  liquid  is  completely  defined 
by  a  numerical  determination  of  a  such  that 

K  -  K  (2-41) 

N  R 

The  value  of  o  for  which  this  condition  is  satisfied  may 

be  defined  as  the  consistent-vibrations  parameter,  a  .  Once 

c 

a  has  been  determined  by  using  trial  values  of  a  in  the 

C  *' 

evaluation  of  Equations  2-33  and  2-40,  the  well  depth  can 
be  computed  by  use  of  2-37.  Expressions  for  the  remaining 
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quantities  which  are  required  in  the  evaluation  of  thermodynamic 
properties  according  to  the  present  model  are  derived  below. 

2.3  Internal  Energy  and  Entropy  of  the  Cuiescent  Liquid 

In  a  system  of  !I  molecules  the  total  energy  of  vibration 
(including  the  “zero-point"  energy)  is 


r 


F 

-4 

J 


r 


j 


(2-42) 


where  E ,  and  F 

3  j 

respectively. 


are  given  by  Equation  2-16  and  Equation  2-19, 
A  series  identity  gives  for  one  mole  of  liquid 


A 


E  ■  RV(e,T) 


(2-43) 


where  R  is  the  universal  gas  constant  and  where 

V(9,T)  =  36  |  |  ♦  (  e9/T-l)  ]  (2-44) 


In  the  present  mode,  fis,  of  course,  computed  by  substituting 
the  consistent  a,  o^,  into  Equation  2-17. 

If  the  liquid  were  composed  of  N  stationary  molecules 
confined  to  their  central  lattice  sites  by  wells  of  depth  u, 
the  internal  energy  of  the  system  would  be  (1/2)  Nu.  This 
result  follows  from  the  fact  that  the  energy  which  belongs 
to  each  molecule  uniquely  is  only  j  u,  the  remaining  half 
figuring  into  the  potential  energies  of  neighboring  molecules. 
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In  other  words,  when  one  adds  up  **«  potential  eaer-ry  of 
U*  syiter,  the  interaction  between  each  pair  should  be  courted 
only  once.  Hence,  the  molar  internal  energy  of  tne  quiescent 

liquid  is 

-  i  J-  ♦  3V(f,7)l  (2-45) 

a*.  » 

The  entropy  of  a  lattice  conposed  of  X  harmonically 
bound  molecules  and  n  vacancies  is  given  by 


S  *  S  ♦  S 
Let 


(2-46) 


where  Sc  and  are  the  configurational  and  thermal  entropies, 
respectively.  According  to  the  Boltzmann  postulate  these 
are  given  by 


S  -  Sc  in 
c  c 


(2-47) 


and 


S  »  k  in  f:t 


(2-43) 


where  W  is  the  number  of  distinguishable  ways  that  the  K 
molecule  nay  occupy  to  N*-n  lattice  sites  and  ;J  is  the  nuiaber 
of  distinguishable  ways  that  the  N  molecules  nay  populate 
the  permissible  quantum  levels  in  a  systea  of  tenperature  T. 
Thus,  is  the  number  of  ways  that  n  identical  pegs  can 
occupy  N+n  holes,  namely. 


*7  .  <w+n>  1 

’  c  Nlnl 


(2-49) 
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For  a  hole  concert rat ice,  f,  ieflaei  by 


(2-s:) 


use  of  the  darinamt  terms  cf  Stirling"*  approximation  give* 
Tor  the  configuration  er.tropy  of  or*  -ole  of  liquid 


S  »  -5  f  m(l-f) 

c  l 


(2-51) 


The  quanta  which  deterring  the  rusher  of  distinguishable 
"energy  arrangements,"  '.'t,  cerrespoed  to  the  total  internal 
energy  is  excess  of  the  grcusc-atate  energy.  In  other  words, 

-■»  cannot  include  pensatatica  involving  depletions  of  the 
ground  states.  Pror  Equation  2-16,  Equation  2-43,  and 
Equation  2-44  it  is  sees  that  the  number  q  of  quanta  which 
contribute  to  the  energy  is  excess  cf  that  cf  the  ground 
states  oust  satisfy  the  relation 

q  k  f  «  ::(3k«)  (e#/T-lj  (2-52) 


from  which 


q/3N  »  g  *  (  e9/T-l  ) 


-1 


(2-53) 


where  g  has  been  introduced  for  sirplicity.  The  number  of 
distinguishable  ways  that  the  q  quanta  can  be  distributed 
among  the  3K  oscillators  is  the  same  as  the  number  of  dis- 
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ti^sitiubie  vtys  that  q  identical  balls  nay  be  pot  is  3?: 
ocxes  wits  so  restriction  os  the  number  of  balls  per  box, 
sareiy, 


-<  «  (3r+i)  1/(321)  Sql  (2-54) 

t 

T-ier,  by  Stirling's  approximation  the  therral  entropy  per 
scle  is 

S  »  3R  ( ( l*c)  in  ( l-*g)  -g  ta  g)  (2-55) 

The  foregoing  equations  define  the  thensodynaric  properties 
of  the  proposed  quiescent  liquid  with  a  hole  concentration  f. 

2.4  Acoustical  Effects 

According  to  Assumption  7,  the  quiescent  sys Uc  is  disturbed 
by  acoustical  waves  with  a  velocity  determined  by  the  ■consistent" 
frequency  of  vibration  and  with  wavelengths  restricted  to 
integral  multiples  of  twice  the  lattice  dissension.  The  average 
tire  required  for  a  siolecule  to  cosrunicate  the  approach 
of  a  disturbance  across  the  distance  2a,  which  separates  its 
nearest  neighbor  on  opposite  sides,  is  1/v,  where  v  is  the 
consistent-vibration  frequency.  The  propagation  velocity 
in  a  consistently  vibrating  liquid  (and  in  a  region  devoid 
of  holes)  is  thus  taken  to  be 

w  =  2av  *  2ak6/h  (2-56) 
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wh*x*  a  is  identical  to  the  nearest-neighber  separation 

in-  where  •  is  obtained  by  substituting  the  consistent  vibration 

parameter,  «  ,  into  2-17. 
c 

The  permissible  wave  lengths  are  given  by 


X  «  n(2a» 
n 


(2-57) 


where  n*  1,  2,  3,....  The  permissible  acoustical  frequencies 
are  therefore  given  by 


r.  "  w/*n  *  *e/hr- 


(2-58) 


Since  it  is  to  be  expected  that  the  acoustical  disturbances 

do  not  extend  over  conains  which  include  vacancies  in  the 

qua  si -lattice,  the  number  of  waves  of  a  given  length  should 

decrease  rapidly  with  increasing  hole  concentration  f.  The 

situation  can  be  approximately  represented  by  a  simple  cubic 

lattice  consisting  of  N  sites  of  which  fN  are  unoccupied. 

1/3 

The  length  of  a  linear  column  of  molecules  is  N  a,  which 


neans  that,  in  the  absence  of  holes,  the  number  of  one-dimensional 

•1/3/2„. 


waves  of  length  X  that  can  exist  alone  the  column  is  N 


However,  if  each  site  has  a  probability  f  of  being  unoccupied, 

the  average  number  of  one  dimensional  waves  that  can  exist 

2n  1/3/ 

along  the  column  is  only  (1-f)  N  /  2n.  Considering  that 
2/3 

there  are  N  columns  for  each  of  three  mutually  orthogonal 
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directions,  the  number  of  one-dimensional  acoustical  oscillators 
of  frequency  in  one  mole  of  liquid  is  (taking  into  account 
the  increase  in  the  number  of  sites  per  mole  due  to  the  vacancy 
fraction  f) 


N  « 
n 


(2-59) 


According  to  Assumption  7,  the  acoustical  contribution 

to  the  internal  energy  is  obtained  by  multiplying  by  the 

energy  of  a  one-dimensional  oscillator  of  frequency  v  averaged 

n 

over  states  with  energies  less  than  k0  and  not  including 

the  ground  state  energy.  Since  the  assured  absorption  of 

waves  of  energy  ke ,  or  greater,  merely  changes  the  form  of 

the  available  kinetic  energy,  the  level  population  of  oscillators 

of  frequency  v  should  not  be  affected  by  the  process  of 
n 

truncation.  Hence,  the  average  energy,  including  the  ground 
state  energy,  is 


where  F  .  is  the  fraction  of  oscillators  of  frequency  v 
n  j  n 

which  are  in  the  jth  level  in  an  ordinary  quantized  distribution 
characteristic  of  a  temperature  T,  E  .  is  the  energy  of  the 
jth  level,  and  J  is  the  level  corresponding  to  an  excitation 
ke.  For  a  one-dimensional  oscillator, 
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(2-61) 


F  *  c 
nj 


jhv  /kTI  -hv^T 
v  'l-e 


2nd 


F  «  -  hv  +  jhv 

nj  2  n  J  n 


The  first  summation  in  the  numerator  of  Equation  2-60 


(2-62) 


is 


E 

3-o 


F  ,E  .  * 
nj  nj 


j-0 


m 

I  -hv  AT\  r-’  -hv  /kT 
4  J  l-e  **  |hvn  )  je  r* 


I  hvn  +  h“,/l 


j-0 


hVn/kT 

e  -1 


(2-63) 


The  second  summation  in  the  numerator  of  Equation  2-60  can 


be  written 


I 

j-J 


FnjEnj 


IhvJl-e‘h'''/kT)£e"^nAT 

j-J 


(2-64) 


The  foregoing  expression  nay  be  simplified  by  writing 
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•  m 

-jhv ^ct  -Jhv^cT  ^  -(j-J)hv 


j*j 


j-j 


•  z« 


-Jhv^T  Y"  -j'h'-./tT 
€ 

j  '*0 

-Jhv 


vkT|i.e-hVTf 


(2-65) 


and 


-jhv^ 


E  -jhv  AT  r-  -)hv  AT  r-* 

je  V  -  )  (W)e  ♦  J  )  e 

j-j  j-j  j-j 


-K'tA' 


-Jhv 


**TZ 


-jCh  At 

ye  * 


j'-0 


♦  Je'Jhv'/tT  £  e-rh“^ 


j'-0 


-Jhv  ^AT 


-h.^t  i 


( .-.-v" 

-jhv  At  I  -hv  A 
+  Je  V  ll-e  ** 


-1 


(2-66) 

Combining  Equation  2-64 .  Equation  2-65,  and  Equation  2-66  gives 


L  V»J 


j*J 


r 

/  \~l  " 

’ 

a  1 

1  ■JhvI/'tT  -Jhv^T 

■r  hv  e  r  thv  e  f 

> 

2  n  n 

- 

j 

f  (2-67) 
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Combining  Equation  2-63  and  Equation  2-67  gives  for  the  numerator 
of  Equation  2-60 


f 


[ 


x  hv  1-e 
2  n\ 


f  F  e  -  t  r  E 

Z— •  nj  nj  (— j  n} 

j*C  :«J 
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The  summation  in  the  denominator  of  Equation  2-60  is 
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and  the  denominator  is 


V 

1  -  L  Fnj 

j-J 


-Jhv  yk: 

1  **€■ 


From  Equation  2-68  and  Equation  2-70,  it  is  found  that 


(2-63) 


(2-69) 


(2-70) 
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(2-71) 


Since  the  Jth  level  corresponds  to  the  absorption  of 
an  amount  of  energy  ke  by  a  ’consistent"  node  of  vibration. 


Jhv  ■  ke 
n 


(2-72) 


and 


J  -  ke, 


e/hv  *  n 
/  n 


(2-73) 


where  use  has  been  made  of  Equation  2-58.  Thus, 


1  .  Ake 
u  »  -  hv  ♦  — 
n  2  n  n 


e‘/n-l  e*-l 


(2-74) 


where 


♦  -  8/T 


According  to  Assumption  8  the  ground-state  energy  of  the 
oscillators  in  normal  coordinate  space  is  to  be  identified 
collectively  with  the  energy  of  the  consistently  vibrating 
molecules.  Hence,  the  average  energy  belonging  uniquely 
to  an  oscillator  of  frequency  is 

u'  *  ~  [ F  -nF  1  (2-76) 

n  n  L  n  1J 
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where 


n 


L,/n-il 


-1 


(2-77) 


The  total  acoustical  energy  per  mole,  U  ,  is  obtained  by 

A 

sunning  over  all  wavelengths: 

n«I  n*l 


where 


H 


(1-f) 


2n 


n  (2n) 


(2-78) 


(2-79) 


A  sort  of  configurational  entropy  is  associated  with 

each  set  of  oscillators  of  frequency  owing  to  the  one- 

dinensional  regions  over  which  waves  are  excluded  due  to 

the  presence  of  vacancies.  The  number  of  intersections  of 

intervals  over  which  waves  could  extend  in  the  absence  of 

holes  is  (1-f)^  N1/3/2n.  The  average  number  of  holes  in 

1/3 

such  an  array  i3  .  Since  the  total  range  excluded  by 

holes  may  be  between  two  wave-supporting  regions,  there  is 
effectively  no  limit  on  the  number  of  vacancies  which  separate 
adjacent  waves.  Hence  the  number  of  distinct  combinations 
in  a  one-dinensional  array  of  N1//3  molecule  is  (in  analogy 


to  Equation  2-54) 

t 


W' 
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(2-80) 
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Considering  that  the  assignment  of  a  given  vacancy  distribution 
2/3 

to  each  of  J!  '  columns  completely  determines  the  vacancy 
distribution  in  the  system  of  N  sites,  the  total  number  of 
distinct  combinations  is 


\1 

A,n 


,2/3 


W' 


A  ,n 


and  the  associated  configurational  entropy  per  nole  of  oscillators 
of  all  permissible  frequencies  turns  out  to  be 


iirr' 

1 

1*7! 

1 1+G  ' 

j tnl 1+G  j-G  InG  1 

1  Z— 

l  n 

11  n|  n  nj 

r.“l 


(2-81) 


where 


G 

n 


(l-f)2n  2ni:n 
2nf  *  f 


(2-82) 


Since,  due  to  reflection  from  holes  or  boundaries,  a 
given  acoustical  wave  may  propagate  to  any  point  in  the  liquid, 
various  waves  of  the  same  frequency  are  indistinguishable 
from  one  another  in  the  sense  that  there  is  no  way  of  deter¬ 
mining  whether  a  wave  corresponding  to  a  given  excited  state 
of  oscillation  is  the  sane  wave  which  was  excited  in  a  previous 
condition  or  a  different  wave  which  has  come  to  occupy  the 
position  of  the  former  one.  Hence  it  is  not  meaningful  to 
speak  of  the  distinguishable  ways  in  which  quanta  of  energy 
may  be  dirtributed  among  the  abstract  oscillators  in  normal 
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coordinate  space  that  define  a  collection  of  waves  of  frequency 
v  .  The  thermal  entropy  associated  with  such  distributions 
is  therefore  zero.  However,  the  entropy  given  by  Equation  2-31 
could  be  regarded  as  a  kind  of  thermal  entropy,  since  it 
relates  not  just  to  the  number  of  distinguishable  hole  dis¬ 
tributions  but  to  the  permissible  patterns  of  motion  for 
acoustical  disturbances  in  the  liquid. 

The  contributions  to  the  total  entropy  denoted  by  Equation  2-51, 
Equation  2-55,  and  Equation  2-31  are  the  only  ones  to  be 
considered  in  a  consistent  application  cf  the  present  rodel. 

An  omission  which  should  be  noted  explicitly  i3  the  part 
of  the  configurational  entropy  associated  with  long-range 
disordering  in  the  distribution  of  quasi-lattice  sites. 

It  is  to  be  hoped  that,  well  away  from  the  critical  point 
at  least,  this  contribution  is  small.  In  the  vicinity  of 
the  critical  point,  any  model  based  on  the  assumption  of 
a  crystalline  structure  can  be  expected  to  fail  (Ref.  11) . 

2,5  Hole  Concentration  and  the  Approach  to  Equilibrium 

According  to  Assumption  9,  the  hole  concentration,  f, 
is  the  only  significant  independent  variable  in  the  determina¬ 
tion  of  the  Helmholtz  free  energy.  Therefore,  at  constant 
temperature  and  density,  the  liquid  should  approach  equilibrium 
by  tending  toward  a  value  of  f  for  which  the  Helmholtz  free  energy 
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F  »  U-TS 


(2-83) 


is  minimum,  where  U  snd  S  are,  respectively,  the  total  internal 
energy  and  entropy  due  to  the  various  contributions.  The 
equilibrium  state  of  the  liquid  is  thus  determined  by  (1)  regard¬ 
ing  each  of  a  set  of  trial  values  of  f  as  an  independent 
variable  (along  with  0  and  T)  and  determining  the  corresponding 
consistent-vibration  frequency  defined  by  Equation  2-41;  (2)  evalu¬ 
ating  the  associated  values  of  U  and  S  by  adding  0^  and  and 

S  ,  and  S  ,  and  Sm ,  respectively;  and  (3)  locating  the  value 
t  c  A 

of  f  for  which  F  is  minimum.  All  of  these  operations  can 
best  be  performed  numerically  with  the  aid  of  a  high-speed 
computer. 


2.6  Summary  of  Expressions  for  the  Equilibrium  State 


The  equilibrium  state  of  a  liquid  according  to  the  present 


model  is  defined  by  an  equilibrium  vacancy  fraction  f#  and 

a  consistent-vibration  parameter  ac|f#j.  For  any  trial  value 

of  f,  a  (f)  is  the  value  of  a  that  satisfies 
c 


K  -K  »  0 

N  R 

where 

\ 

and 

■  Z  h  [*rSk«lYi)  *  \"14|‘r(,i)] 
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(2-41) 

(2-40A) 

(2-33) 
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(2-34) 


(2-35) 


(2-36) 


The  polar  well  depth 


*•  W  A. 


given  by  substituting  a  *  a  (f)  into  the  relations 
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du  (2-3S) 
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The  lower  Unit  in  the  integrations  is  taken  to  be 


L  (Ti)  ■  i[n  */iJ"26  ]  ( 

The  nunber  of  neighbors  per  shell,  i. ,  and  the  neighbor  shell 

th 

radius,  a. ,  for  the  i  *  shell  are  proportional  to  (1-f)  and 
x 

(1-f)1^3,  respectively,  and  are  given  for  sone  specific  types 
of  quasi-lattice  structure  in  Table  1. 

The  internal  energy  for  a  trial  value  of  f  is 


U  »  U  +  U 
C  A 


where 


with 


and 


V( 0 ,T)  *  30 


and  U  is  given  by 
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(2-79) 


and 
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(2-77) 


F  «  (  e*/n-l 

n 


-1 


The  absolute  entropy  for  a  trial  value  of  f  is 


S  «  S  ♦  S  +  S 
c  t  A 


where 


S  *  -R 
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with 
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The  Helmholtz  free  energy  for  a  trial  value  of  f  is 

F  «  U  -  TS 

The  equilibrium  state  of  the  liquid  is  found  by  locating 
a  value  of  f  such  that 


il 

af 


P  ,7 


(2-85) 

(2-51) 

(2-55) 

(2-53) 

(2-31) 

(2-82) 

(2-83) 

(2-36) 
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If  it  is  assumed  that  the  present  model  provides  accurate 
values  of  U  and  S,  then,  for  a  given  density  p  and  temperature 
T,  the  pressure  is  obtainable  from  the  relation 


Therefore,  insofar  as  it  is  valid,  the  model  offers  a  complete 
solution  to  the  problem  of  state  with  p  and  T  as  independent 


variables. 


3.  THE  COMPUTER  PROGRAM 


3.1  Procedure 

A  conputer  program  (designated  T-74)  has  been  written 
in  the  FORTRAN  language  for  use  on  the  IBM  7090  conputer. 

Tne  operations  indicated  in  Section  2.6  are  executed  on  tne 
basis  cf  given  temperature,  density,  and  assured  quasi-lattice 
structure.  Up  to  ten  neighbor  snelis  ray  be  considered.  For 
each  of  three  trial  values  of  f,  provided  as  input,  the  procedure 
converges  on  a  consistent-vibration  parameter,  a  (f) ,  by  an 
iterative  process  equivalent  to  Newton’s  rsethod,  which  begins 
with  the  consideration  of  two  trial  values  of  a,  also  provided 
as  input.  Subsequently,  the  procedure  converges  on  an  equili¬ 
bria r.  value  of  f  by  re-use  cf  Nevcon’s  method  to  locate  the 

intercept  of  OF/Jfj^  .  A  linear  expression  for  the  latter 

*■  #3 

quantity  at  each  value  of  f  considered  is  obtained  by  using 

th*»  derivative  of  the  quadratic  function  defined  by  the  values 

of  F  at  that  value  of  f  and  at  the  two  values  considered  irrcdi- 

ately  prior  to  it.  In  order  tc  ootain  five- place  accuracy 

in  %  ,  about  six  trial  values  of  a  rust  be  considered  for  eacn 
c 

trial  value  of  f.  In  order  to  obtain  four-place  accuracy  ir. 

f  ,  about  1")  trial  values  of  f  oust  re  considered.  (For  sere 
e 

reason,  the  number  of  iterations  required  in  each  case  is  rela¬ 
tively  ir sensitive  to  the  estirates  provided  as  input.)  The 
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execution  of  all  of  the  foregoing  operation?  for  a  given  state 
defined  by  T  and  p  (in  the  cases  considered  thus  far,  at  least) 
requires  about  3  minutes  of  computer  time. 

3.2  Input  and  Output 

The  input  data  for  a  single  problem,  (corresponding  to  a 
single  state  of  the  liquid)  requires  from  7  to  16  cards,  depend 
ing  on  the  number  of  neighbor  shells  considered.  If  n  is  the 
number  of  neighbor  shells,  then  the  information  contained  on 
the  respective  input  cards  is: 


11)  Twc  acbreviations  designating  the  type  of 
liquid  and  assured  structure  (used  in  the  identi¬ 
fication  of  input) i  the  number  of  neighbor  shells, 
n;  the  temperature?  the  density  (used  only  to 
specify  what  density  is  used  in  determination  of 
the  shell  radii  for  the  perfect  lattice) ?  the 
molecular  mass  ir.  amu. 


(2)  The  Lennard- Jones  coefficients  uc  and  v 
(in  units  of  ergs  x  10“®  and  ergs  x  10"®, 
respectively) . 


13  through  3+(n-l)J  The  radius  and  number  of 
molecules  per  shell  for  each  neighbor  shell  of 
a  perfect  lattice  of  the  assumed  type  (from  this 
information  the  corresponding  values  in  a  lattice 
with  a  hole  concentration  f  are  computed  by 
multiplication  by  the  factors  (1-f)^^  and  (1-f)  , 
respectively t  thus,  p  is  implicit  in  the 
■perfect-lattice"  input) . 


(3+r.J 


Two  trial  values  of  a 


(3+n+l]  Three  trial  values  of  f. 


l3+n+2]  The  mesh  spacing  (£u)  to  be  used  in  the 
integrations  and  an  optional  lower  limit  of 
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••m/ntfumt-i- 


If 


integration  which  can  be  used  in  the  case  of  distant 
neighoor  shells  to  reduce  the  amount  of  machine 
time  required.  The  upper  limit  of  integration 
in  all  cases  is  automatically  chosen  to  give  equal 
ranges  of  integration  on  either  side  of  the  maxi¬ 
mum  of  the  integrand  associated  with  k_r  J  y ^  j  • 

(3+n+3)  An  option  number  which  permits  tne  termi¬ 
nation  of  the  convergence  process  and  which  also 
permits  the  inclusion  or  exclusion  of  acoustical 
terms  (the  properties  of  the  "quiescent"  liquid 
are  of  some  interest) ;  two  integers  which  specify 
an  upper  limit  to  the  number  cf  trial  values 
of  a  and  f  to  be  considered  in  a  given  problem; 
two  (fractional)  numbers  which  define  the  con¬ 
vergence  criteria  and  hence  the  accuracy  of  oc 

and  f  . 
e 

The  output  data  include  all  of  the  input  data  and,  for 
each  trial  v. lue  of  f,  the  following  in  the  order  considered: 

1.  The  value  of  f  considered 

2.  The  "consistent"  values  of  y^  for  each  neighbor 
shell 

3.  The  contribution  to  the  "consistent"  stiffness 
coefficient  due  to  each  neighbor  shell 

4.  The  contributions  tc  the  effective  molar  well 

depth  (designated  as  U  *=  1/2  W)  due  to  each 
neighbor  shell  g 

5.  The  total  effective  molar  well  depth  and  stiff¬ 
ness  coefficient 

6.  The  consistent  values  of  u,  which  is  a  ,  and  9 

7.  The  total  internal  energy,  total  entropy,  m'1 
the  Helmholtz  free  energy 

S.  The  value  of 
sidered  r 


JL 

PT 


JF 

if 


for  tae  val’ie 


of  f  '~on- 
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Each  value  of  f  corresponds  to  a  single  page  of  output j 
in  a  normal  execution  the  last  page  of  output  provides  the 
above  data  for  the  equilibrium  state  of  the  liquid. 
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4.  ATTLICATIOX  TO  LIQUID  HYDROGEN 


4 . 1  Liquid  Structure  and  Intermolecular  Potential 

At  first  glance,  the  application  of  the  consistent-vibrations 
model  to  liquid  hydrogen  night  seen  inauspicious  due  to  the 
fact  that  hydrogen  is  diatonic.  However,  infrared  spectroscopy 
(Ref.  12)  reveals  that  the  rotation  o*  the  molecule  in  the 
liquid  (and  even  in  the  solid)  is  practically  unhindered,  and 
analysis  of  virial-coef ficient  data  on  hydrogen  gas  shew  that, 
when  both  the  para-  and  ortho-hydrogen  modifications  are  taken 
into  account,  the  interaction  between  molecules  is  approxi¬ 
mately  satisfied  by  a  Lennard-Jones  potential  (Ref.  13) .  Further¬ 
more,  in  either  modification  the  first-accessible  excited 
rotational  level  of  the  molecule  lies  approximately  0.03  ev 
above  the  ground  state.  Hence,  even  at  temperatures  as  high 
as  30*  X,  the  fraction  of  molecules  which  are  in  the  first 
excited  state  of  either  modification  is  only  of  the  order 
exp(-.03/kT)~  10  5.  Thus,  contributions  to  the  internal  energy 
due  to  excited  rotational  states  may  conscientiously  be  neglected. 

The  Lennard-Jones  parameters  which  have  been  used  in  all 
calculations  to  date  are  those  cited  by  Knapp  and  Beenakker 
(Ref.  13).  These  are  (in  the  notation  and  units  selected  *or 
the  present  application) 

u  =  0.001257  A6 ( ergs  x  10~q  )  (4-1) 
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and 


C'tfilj 


•  I 

A|«rg* 


x  10 


-8) 


(4-2) 


These  parameter*  are  obtained  from  the  properties  of  hydrogen 
gas  and  are  not  sufficiently  accurate  to  justify  a  distinction 
between  the  ortho-  and  para-hydrogen  modi ficat ion*  (although 
the  relative  difference  between  the  coefficients  of  the  two 
modifications  can  be  specified,  as  is  done  in  Reference  13) . 

The  structure  of  the  quasi-lattice  in  liquid  hydrogen 
is  assumed  to  be  cubic  c.ose-packed.  In  all  applications  made 
thus  far,  the  first  four  neighbor  shells  have  been  considered . 
The  possibility  that  the  type  of  lattice  structure  corresponding 
to  the  lowest  Helmholtz  free  energy  ray  change  from  cubic  close- 
packed  to  some  other  form  at  the  higher  temperatures  has  yet 
to  be  investigated. 

4.2  Boult.  .^.Ccpg.rlao,  .1th  Ey rlont.1 

The  procedure  has  been  applied  to  the  case  of  liquid  para- 
hydrogen  at  points  on  the  P-V-T  surface  ranging  from  20 *R 
to  28*K  along  the  saturation-point  line.  Typical  results 
are  exemplified  by  the  output  for  the  state  defined  by 
T  ■  20.26C*  X ,  p  *  0.0707?  gm/cra3,  which  is  shown  in  Table  2. 

The  data  refer  to  the  equilibrium  state  of  the  liquid.  In 
the  case  shown,  the  contribution  of  acoustical  terms  was  not 
considered. 
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Table  2 


I 


T-74  OUTPUT  FOR  THE  STATE  OF  LIQUID  PARA-HYDROGEN 
(t  *  20.268*K,  p  *0.07077  gm/cm3) 


i 

yi 

KP.(i)*j 

ergs/cn2] 

1/ 2T 7  < i )  *  ( j/gm) 

1 

6.6489 

322.049 

-453.401 

2 

9.4029 

'.914 

-39.372 

3 

11.5163 

—6 

i  .622 

-45.108 

4 

13.2978 

-] 

..042 

-9.420 

Ho.  of  o*s  considered  *  6 
(Ho.  of  f's  considered)  *  8 

(Equilibrium)  Hole  Concentration  jf^j  *  0.0S717 
Consistent  o  jac)  *  1.7392  A-1 
(Consistent)  0  *  73.388*K 

(Effective)  Hell  Depth  (1/2H)  *  -546.301  j/gm 
Internal  Energy  juoj  of  Quiescent  Liquid  =  -304.836  j/grr 
Entropy  °~  ~uiesc©nt  Liquid  *  2.9284  j/cr-*K 

Helmholtz  Free  Energy  of  Quiescent  Liquid  *  -364.189  j/gm 


RT 


*0.00668 


•Contributions  to  KR  and  1/2  W  from  the  it*1  neighbor 
shell. 


1€ 


All  of  the  problems  referred  to  herein  were  run  before 
the  computer  program  had  been  modified  to  include  the  contri¬ 
bution  of  U  to  the  internal  energy  or  the  contribution  of  S 
A  A 

to  the  entropy.  The  subsequent  evaluation  of  these  quantities 

by  hand  has  yielded  values  for  the  total  energy  and  entropy 

which,  in  terms  of  the  model,  are  valid  for  the  equilibrium 

hole  concentration  found  by  the  computer  program  but  which 

do  not  necessarily  correspond  to  the  vacancy  fraction  that 

would  have  been  determined  had  U.  and  £.  been  included  in  the 

A  A 

minimization  of  F.  Nevertheless,  since  the  hand  calculations 
show  that  U  represent  only  about  20%  of  the  absolute  value 

A 

of  U  amd  that  varies  slowly  with  f,  it  is  inferred  that 

the  value  of  the  equilibrium  vacancy  fraction  will  not  be  greatly 

disturbed  by  the  inclusion  of  such  terms. 

The  most  significant  results  of  the  computer  calculations 

and  the  hand  calculations  of  U.  amd  Sm  are  shown  in  Tad>le  3. 

A  A 

It  is  to  be  noted  that  f  increases  monotonically  (and  almost 
linearly)  as  one  moves  toward  the  critical  point  (33*K) .  This 
is  just  what  would  be  expected,  since  all  quasi-crystalline 
ordering  is  supposed  to  vanish  in  the  neighborhood  of  the  critical 
point.  The  quantities  e.-  which  vary  only  slightly,  are  of 
immediate  value,  since  they  determine  the  liquid- lattice  quantum 
levels  that  may  be  excited  by  ccld-neutron  scattering.  The 
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internal  energy  and  entropy  of  the  quiescent  liquid  are  seen 
to  depend  rather  strongly  on  the  state.  As  remarked  previously, 
the  acoustical  contribution  to  the  total  entropy  appears  to 
be  insensitive  to  the  state,  and  hence  to  fft. 

The  applicability  of  the  proposed  model  can  presumably 
be  tested  adequately  by  comparing  the  predicted  values  of  the 
total  internal  energy  and  entropy  with  empirical  results. 

Problems  involving  variations  of  the  density  at  constant  tern" 
perature  and  vice  versa  have  net  yet  been  run,  so  neither  pre~ 
dictions  of  pressure  nor  of  specific  heat  at  constant  volume 
are  yet  available. 

Figure  2  shows  the  predicted  value  of  the  total  internal 
energy  |uo  ♦  O^J  compared  with  values  published  by  IJ3S  (Ref.  14). 
Figure  3  shows  the  predicted  entropy  compared  wit!  the  empirical 
results  from  the  sane  reference. 
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Figure  2  Calculated  and  Empirical  Internal  Energy  of 
Liquid  Hydrogen  Along  the  Saturation  Line 
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Figure  3  Calculated  and  Empirical  Entropy  of  Liquid 
Hydrogen  Along  the  Saturation  Line 
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5.  RESULTS  AND  CONCLUSIONS 

The  predicted  internal  energies  and  entropies  lie  within 
6%  of  the  empirical  results  at  points  -on  the  saturation  line 
below  about  25*  (Fig.  3) .  At  points  nearer  the  critical  point, 
the  predicted  results,  especially  those  for  the  entropy,  give 
less  satisfactory  agreement.  However,  this  is  to  be  expected 
since  the  concept  of  a  close-packed  quasi-lattice  tends  to 
lose  its  meaning  when  the  equilibrium  vacancy  fractions  approaches 
20t  • 

r.?hen  it  is  considered  that  the  agreement,  in  the  case 

of  the  internal  energy,  at  least,  would  almost  certainly  improve 

if  more  than  four  neighbor  shells  were  considered  [ judging 

from  the  relative  contributions  to  1/2  W  and  K  | ,  the  results 

obtained  thus  far  constitute  highly  encouraging  evidence  *cr 

the  validity  of  the  proposed  model.  It  is  therefore  tentatively 

concluded  that  the  model  anticipawas  the  most  significant  general 

features  of  notion  in  the  liquid  and  t’:at  the  state-dependent 

nicroscopic  data  provided  by  the  model,  such  as  0  and  the  neighbor 

wave  functions  determined  by  o  ,  may  be  implemented  in  calcula- 

c 

tions  of  quantities  such  as  cross  sections  and  recombination 
rates  which  are  dependent  on  the  details  of  the  internal  dynamics 
and  microscopic  structure. 
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